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SUMMARY 


Experimental results are presented for transonic flow past cone- 
cylinder, axially symmetric bodies. The drag coefficient and surface 
Mach number are studied as the free-stream Mach number is varied and, 
wherever possible, the experimental results are compared with theoret- 
ical predictions . Interferometric results for several typical flow 
configurations are shown and an example of shock-free supersonic-to- 
subsonic compression is experimentally demonstrated. 

The theoretical problem of transonic flow past finite cones is dis- 
cussed briefly and an approximate solution of the axially symmetric 
transonic equations, valid for a semi-infinite cone, is presented. 


INTRODUCTION 


Transonic flow past certain two-dimensional bodies has been the sub- 
ject of several recent papers and the phenomena are well understood. The 
theoretical results of Cole (ref. 1 ), Guderley and Yoshihara (ref. 2), 
Vincent! and Wagoner (ref. 3 ) , and others apply to two-dimensional sym- 
metrical double-wedge airfoils. The experimental results of Bryson 
(ref. A) and Griffith (ref. 5) substantiate the theoretical work in a 
very satisfactory manner. More recently, Vincenti and Wagoner (ref. 6 ) 
and Guderley and Yoshihara (ref. 7 ) have analyzed the transonic flow 
past two-dimensional unsymmetrical sections, that is, lifting double- 
wedge airfoils. Current experiments on lifting double -wedge airfoils 
(ref. 8 ) at the Guggenheim Aeronautical Laboratory of the California 
Institute of Technology indicate that agreement between theoretical and 
experimental results will again be obtained. 

Two-dimensional and axially symmetric transonic flows are of con- 
siderable theoretical and practical interest since these two specialized 
problems are limiting cases of the more complex problem of the flow about 
an arbitrary three-dimensional body. 

The study of axially symmetric transonic flow is not so complete as 
that of two-dimensional flow. In recent years several papers, notably 
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those of Von Karman (ref. 9 ) and Oswatitsch and Berndt (ref. 10), have 
studied the similarity rules of axially symmetric transonic flow. Also, 
Yoshihara (ref. 11) has calculated the flow about a finite cone at a 
free -stream Mach number of 1 by a relaxation technique and has obtained 
some experimental verification of the theoretical result. The hodograph 
problem for general transonic flow past finite cones is discussed in 
reference 8 . However, theoretical solutions or experimental results for 
the complete transonic regime are not, at present, available. The pres- 
ent paper presents the results of an experimental investigation of the 
transonic flow past cone -cylinder bodies. A conical tip followed by a 
cylindrical afterbody was chosen as the experimental model for two pri- 
mary reasons: (l) The relatively simple geometry of a cone-cylinder 

body may simplify the theoretical problem, and ( 2 ) viscous effects are 
minimized; that is, the boundary layer on the cone surface is in a 
region of decreasing or constant pressure so that the presence of the 
boundary layer will not greatly alter the shape of the body forward of 
the cone shoulder. 

Theoretical results for the supersonic flow past a cone were first 
presented in 1929 by Busemann (ref. 12). Busemann's solution postulates 
a semi-infinite cone and assumes that the flow is conical; that is, 
along rays through the apex of the semi -infinite cone, the flow parame- 
ters such as pressure and velocity are constant. The solution is found 
by a geometrical construction in the hodograph plane and it is readily 
apparent that a conical solution exists only so long as a shock wave is 
attached to the cone apex and, therefore, the free -stream Mach number 
is supersonic. It is interesting to note that Busemann's solution pre- 
dicts smooth shock-free compression from supersonic to subsonic flow 
for particular combinations of cone angle and free-stream Mach number. 

The conical solution also shows that for a given cone angle and free- 
stream Mach number Mo, the surface Mach number is always less than the 
Mach number immediately behind the conical shock wave; as Moo decreases, 
the surface Mach number decreases and eventually passes from supersonic 
to subsonic values. As was mentioned by Busemann, the conical solution 
for a semi -infinite cone is completely valid for a finite cone so long 
as the flow is everywhere supersonic, but when the surface Mach number 
is less than sonic the perturbation due to the corner or shoulder of the 
finite cone is propagated forward through the subsonic portion of the 
field destroying the conicity of the flow. Thus, the Busemann solution 
is completely valid for a finite cone only so long as Moo is large 
enough so that the surface Mach number is greater than sonic. 

Taylor and Maccoll (ref. 13 ) in 1933 presented the results of a 
numerical integration of the axially symmetric equations of motion for 
conical flow about semi -inf inite cones and also presented experimental 
verification of their theoretical results. Further experimental veri- 
fication by Maccoll (ref. 14) was published in 1937* Both of the above 
papers noted that deviations of the experiments from the theoretical 
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predictions, notably In the shape of the shock wave, are apparent when 
the surface Mach number is subsonic . 

As was mentioned previously, Yoshihara (ref. 11 ) has computed, by 
relaxation methods, the flow about a cone cylinder at Moo = 1.00 and 
has experimentally verified the calculation. However, theoretical solu- 
tions do not exist for the complete transonic regime. Solutions have 
not been developed for the flow past a finite cone when Mm is subsonic 
or when Moo is between sonic and the value of Mm at which Busemann's 
conical solution becomes valid. Drougge (ref. 15) has computed the flow 
field between a detached shock wave and a finite cone by relaxation 
methods j however, the position and shape of the detached shock wave were 
determined initially from schlieren photographs . 

The experimental results reported in this paper cover several inter- 
esting features of the transonic flow about finite cones. The devia- 
tions of the surface Mach number from the values predicted by conical 
theory are examined for values of Moo such that the flow field is tran- 
sonic in nature . The behavior of the surface Mach number for subsonic 
values of M*,, and as Moo approaches sonic from subsonic values, is 
examined in some detail so that an extrapolation to Moo = 1.00 may be 
made . The above surface Mach number data lead naturally to the evalua- 
tion of the drag coefficient, and experimental values of the drag coef- 
ficient in the transonic regime are presented. 

The physical location of the sonic line in a meridional plane of 
the flow about a finite cone is of considerable interest for a theoreti- 
cal study of the problem of axially symmetric flow. With this fact in 
mind, an interferometric analysis was made at several typical values 
of Mm so as to determine the local Mach number fields about finite 
cones. Several examples of supersonic to subsonic shock -free compres- 
sion are experimentally demonstrated. 

Experimental values of the shock-wave angle at the cone tip, par- 
ticularly at values of Moo where the flow field between the shock wave 
and the cone surface is transonic or subsonic in nature, are presented, 
and a comparison with the values from conical, theory is shown. 

The conical solution for flow about a semi -inf ini te cone demonstrates 
that a conical solution does not exist if, for a given cone angle, 

Mm decreases below a certain minimum M»- This minimum M^ is defined 
to be the M» for which shock-wave detachment occurs for a semi-infinite 
cone. Whether or not the shock-wave -detachment Mach number for a finite 
cone can be determined from conical theory is of considerable theoretical 
interest. Experimental values of the detachment distance of a shock wave 
from a finite cone tip, the distance obviously being zero at attachment, 
have been collected from several sources and the results analyzed in this 
report . 


k 


NACA TN 3213 


The transonic equations of motion and boundary conditions as derived 
by Von Karman (ref. 9 ) for socially symmetric flow require several assump- 
tions as to the relative magnitude of various terms in the exact equations 
of motion and the related boundary conditions . To demonstrate that the 
transonic equations retain the important features of the exact equations, 
an approximate solution of the problem of conical flow about a semi- 
infinite cone has been developed employing the transonic equations and 
boundary conditions. A comparison of the exact Busemann solution and the 
approximate transonic solution is presented in the report. 

The author of this report wishes to express his appreciation for 
their helpful advice and criticism to Drs. H. W. Liepmann, J. D. Cole, 
and A. Roshko of the California Institute of Technology. The investi- 
gation was conducted under the sponsorship and with the financial assis- 
tance of the National Advisory Committee for Aeronautics . 
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EXPERIMENTAL EQUIPMENT 
Wind Tunnel 


The transonic wind tunnel at GALCIT is a continuous -flow wind tunnel. 
For supersonic testing, the test-section Mach number may be continuously 
varied over a wide range by altering the shape of one flexible wall. The 
test-section Mach number is varied by changing the area of a sonic throat 
downstream of the test section when subsonic tests are being performed. 

The test-section width is k inches and the height is 9 inches. The design 
of the flexible test-section wall is discussed in reference l6. 


Interferometer 

The interferometer used in the present investigation is of the Mach- 
Zehnder type . Both light paths of the interferometer are passed through 
the wind-tunnel test section, one beam passing through the flow region 
under analysis and the other beam passing through the undisturbed flow in 
the test section upstream of the model. The fringe shifts due to density 
variations in the boundary layer are eliminated since both beams pass 
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through the boundary layer on the test -section walls and are affected 
almost equally. The growth of the boundary layer between the two beams 
is not compensated by the above arrangement, but the effect is of sec- 
ondary importance. A detailed description of the GALCIT interferometer 
and a very complete bibliography on interferometer construction is given 
in reference 17. 


Models 

The models were conical-tipped brass cylinders of 20°, 25°, and 3O 0 
semiangle. The base diameters were between O.3O and 0.50 inch. Thus 
the Reynolds numbers for the tests, with the base diameter of the models 
as the reference dimension, varied from 55,000 to 143,000. The tips of 
the cones were made as sharp as possible and the maximum tip diameter of 
the dullest of the models was approximately 1/3 percent of its base 
diameter. Also, the models were black-nickel plated to improve the photo- 
graphic definition. The angle of attack and angle of yaw were adjusted 
to zero by equalizing the pressure on the cone surface at four annular 
points . 


EXPERIMENTAL TECHNIQUES 
Interferometry 


An experimental investigation of flow phenomena is facilitated by 
the employment of an interferometer to determine the density fields in 
gaseous (or liquid) flows. The interferometer technique possesses the 
obvious advantage of eliminating the need for placing any type of probe 
into a flow region where the presence of the probe may completely alter 
the undisturbed flow field. A disadvantage is also present, however, 
since the values of density are not immediately available as the test 
is in progress . A more serious disadvantage is the fact that the inter- 
ferometer integrates the density values on its light paths (see appen- 
dix A) and, thus, the measurement of density is not localized but is 
influenced by inhomogeneities in the flow which may be well removed 
from the points of interest . 

Two general types of flow are amenable to interferometric analysis, 
namely , two-dimensional and axially symmetric flow. This paper is con- 
cerned solely with axially symmetric flow analysis. A discussion of the 
method employed to reduce the finite-fringe interf erograms , such as fig- 
ure 1, to density distributions may be found in appendix A. The method 
is essentially that outlined in reference 18. An excellent discussion 
of axially symmetric data reduction is given in reference 19 where sev- 
eral references to earlier papers in the field will be found. 
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These earlier papers are mainly concerned with evaluating the 
interferometer-data reduction techniques for axially symmetric flow by 
investigating the flow about cone cylinders at Mach numbers and cone 
angles where the Busemann conical solution was known to be valid. Ref- 
erence 20 presents some experimental results in the same general flow 
regimes as are investigated in this report. 


Sonic -Line Location by Wave Reflection 

The location of the sonic line in a meridional plane of an axially 
symmetric transonic flow can be determined experimentally by at least 
three distinct methods. The first method is by static-pressure meas- 
urements, the second is by interferometric analysis, and the third is 
that of Mach wave reflection from the sonic line . It should be noted 
that at the point of reflection the Mach wave will be perpendicular to 
the streamline direction through the sonic line. 

To locate the sonic line within the flow about a cone, a small- 
diameter probe was placed in the free -stream flow outside of the cone 
shock wave. The probe was in a position such that the probe shock wave 
pierced the cone shock wave and entered the flow field about the cone. 

The shock waves formed by the probe closely approximate Mach waves at 
large distances from the probe and a typical wave reflection is shown 
in figure 2. The perturbations in the flow about the cone caused by 
the waves do not appreciably affect the position of the sonic line as 
shown by figure 3 where a comparison is made of the location of the 
sonic line as found by interferometric analysis and by the wave -reflection 
method. The probe method is much more convenient than the pressure - 
measurement or interferometric method since the phenomenon may be observed 
with a schlieren system, so that the result is obtained visually. 


Pressure Measurements 

The pressure measurements in this investigation were made either on 
a micromanometer (accuracy of ±0.01 millimeter of mercury) or on a nomo- 
graph Mach meter (ref. 21). 


EXPERIMENTAL RESULTS AND DISCUSSION 
General Flow Characteristics 


An analysis of the flow of a compressible fluid about an axially 
symmetric finite cone, that is, a cone cylinder, indicates that five 
distinct regimes of flow are possible. These regimes are given below. 
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Regime I .- Regime I is subsonic flow at infinity with a region of 
locally supersonic flow downstream of the cone shoulder. A schlieren 
photograph of this type of flow is shown in figure 4. It should be noted 
that an extremely weak shock wave originates at the cone shoulder and 
terminates at the downstream "normal" shock wave. The forked appearance 
of the base of the terminating "normal" shock wave is an illusion caused 
by the axial symmetry of the flow. The light rays which pass near the 
surface of the body in the region of the rearward branch of the "fork" 
also pass through the outer portion of the shock and the spurious rear- 
ward branch is caused by the light-ray deflections in the outer portion 
of the shock wave. A meridional section of the shock wave actually 
includes only the front branch of the fork. 

Regime II .- Regime II is supersonic flow at infinity with a detached 
shock wave and subsonic flow between the shock wave and the cone. Fig- 
ure 1 is a finite-fringe interferogram of this type of flow. 

Regime III .- Regime III is supersonic flow at infinity with an 
attached curved shock wave and subsonic flow between the initial portion 
of the shock wave and the cone. A schlieren photograph of this flow is 
shown in figure 5( a )* Taylor and Maccoll's original paper on conical 
flow (ref. 13) includes a schlieren photograph of an attached curved 
shock. 

Regime IV .- Regime IV is supersonic flow at infinity with an 
attached shock wave and mixed supersonic and subsonic flow between the 
shock wave and the cone. A schlieren photograph of this flow is shown 
in figure 5(t>). 

Regime V .- Regime V is supersonic flow at infinity with completely 
supersonic flow between the attached shock wave and the cone surface. 

The Busemann solution applies in this regime and has been verified experi- 
mentally in references 13, 14, 19, and 22. 


Local Mach Number Contours 

The local Mach number contours in a meridional plane for the flow 
about a 25° semiangle cone are shown in figure 6 for flow regimes II, III, 
and IV. The local Mach number contours for a 30° semiangle cone in 
regime II are shown in figure 7 • These data were obtained by interfero- 
metric analysis as discussed in appendix A. 

The local Mach number contours should be normal to the cone surface 
since the cone surface is a flat boundary and any pressure gradient at 
the surface must be parallel to the flat surface. However, near the 
shoulder of the cone cylinder the surface is curved by the effect of the 
corner expansion on the boundary layer , and thus the local Mach number 
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contours are not quite perpendicular to the cone surface immediately 
ahead of the shoulder. 


Sonic -Line Location 

The location of the sonic line can be determined by interferometric 
analysis, but a more useful method, in the present investigation, was the 
wave -reflection method. The location of the sonic line in regimes II, 
III, and IV is discussed below. 

Regime II .- Examples of the sonic-line location in regime II are 
shown in figures 6(a) and 6(b). The sonic line originates at the cone 
shoulder and terminates on the detached shock wave. In figure 7 it may 
be seen that a region of supersonic to subsonic compression exists on the 
outer portion of the sonic line. The sonic line actually originates 
slightly upstream of the cone shoulder. This effect is due to the 
rounding of the cone shoulder by the surface boundary layer. 

Regime III .- Figure 6(c) illustrates the case of the flow with a 
nearly attached curved shock wave. Again, a small region of supersonic 
to subsonic compression is present on the outer portion of the sonic 
line. The free-stream Mach number is slightly less than the detach- 
ment Moo predicted by the exact conical theory. The question of experi- 
mental detachment Mach number is discussed below. 

Regime IV .- Several examples of the sonic -line location in regime IV 
are shown in figures 3> 6(d), and 8. Figure 3 shows the location as 
determined by interferometric means and as determined by wave reflection. 
The sonic line again originates at the corner and now terminates at the 
cone tip and not on the shock wave as in regimes II and III. A shock- 
free supersonic -to-subsonic compression occurs on the forward portion of 
the sonic line . The location of the sonic line for a 20° semiangle cone 
is shown in figure 8. The agreement between the theoretical and experi- 
mental location is satisfactory near the tip of the cone and for some 
distance downstream of the tip.l 

The question of smooth shock-free supersonic -to-subsonic compres- 
sion has been the subject of much discussion in recent years. The above 
experimental results demonstrate that such a flow is possible. However, 
the smooth compression is not of primary importance, but rather the con- 
ditions under which it occurs. These conditions are that the sonic sur- 
face bounds a zone of subsonic flow completely enclosed by a region of 
supersonic flow and a solid surface . 

1 The existence of this type of flow is indicated by the experimental 
results of Taylor and Maccoll (refs. 13 and 14) and was also discussed by 
Tsien (ref. 23)- 
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As an example of non-shock-free supersonic -to -subsonic compression, 
consider the flow past a two-dimensional airfoil at high subsonic speeds . 
The local supersonic zone on the airfoil is terminated by a shock wave 
and smooth compression through sonic velocity does not occur. In the 
two-dimensional case, however, the supersonic zone is bounded by a sub- 
sonic region and a solid surface . This is the opposite arrangement to 
that in the flow about a cone, in regime IV, where shock -free supersonic - 
to-subsonic compression does occur. 

The above considerations illustrate that the existence (or stability) 
of shock-free compression through sonic velocity may not be a local phe- 
nomenon but may depend on the arrangement of the complete flow field. 


Shock -Wave Angle 

The angle of the attached shock wave at the nose of the cone was 
determined for a 20° and a 25° semiangle cone in flow regimes III, IV, 
and V. The values are shown in figure 9. Similar experiments are 
reported in references 13 and lA. Reference 22 presents data for one 
cone angle at one Mach number in regime III and one Mach number in 
regime IV. The agreement between the exact theory and the experimental 
values at the cone tip is very good even in regimes III and IV where the 
exact theory is not applicable for the complete finite cone. 


Surface Mach Number Distribution 

The distribution of the surface Mach number Mg on a 25° semiangle 
cone for various values of M*, is shown in figure 10 and that on a 20° 
semiangle cone, in figure 11. Several characteristics of these distri- 
butions are of particular interest . 

(1) Surface Mach number near the shoulder deviates from the Busemann 
conical values as soon as M s = 1.00 is attained. Surface Mach number 
near the cone tip agrees quite well with the conical values until the 
theoretical detachment Moo occurs. At the corner M s should, except 
for boundary-layer effects, always be sonic if M s forward of the 
shoulder is subsonic. 

(2) As Mo, approaches 1.00 from the subsonic or from the super- 
sonic regimes. Mg at a particular chordwise station approaches a con- 
stant value. This behavior implies that 


dM„ 


= 0 


Moo=1.00 
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The same behavior of Mg on two-di m ensional sections was noted in 
reference 2b, and thus the concept of stationary values of M s at 
Mo, = 1.00 is established for two-dimensional and axially symmetric flow. 
Since these two cases represent limiting cases of the flow about general 
three-dimensional bodies, the stationary Mg concept can probably be 
applied quite generally if suitable care is taken in choosing the range 
of Moo about Moo = 1.00 in which the so-called "Mg freeze" is 
applicable . 

(3) As Moo progresses from a subsonic value through Moo = 1-00 
and on to a value in regime V, the M s at a particular chordwise sta- 
tion probably varies quite smoothly with no abrupt variations, even at 
attachment of the shock wave, except for a region quite near the tip 
where large variations may occur when the shock wave attaches. 


Drag Coefficients 

The drag coefficients for the 20° and 25° semiangle cones are shown 
in figure 12. The values at M^ = 1.00 were determined by extrapolating 

the M s data in figures 10 and 11 to M^, = 1.00. 


Using the concept of stationary values of M s at M^ = 1.00, the 
drag-curve slope at Mo, = 1.00 becomes (see ref. b and appendix B) 


dC D 




M^l.00 


b 

7+1 


2 

7 + 1 


( 1 ) 


where C^* is the drag coefficient at M^ = 1.00. The first term 

4/(7 + 1) of the drag-curve slope is derived from the first-order term 
of the pressure coefficient while the term 0/(7 + l)J C D * represents 

the contribution of the second-order terms. The magnitude of the 

second-order term { 2/(7 + ljj C^* is shown by the difference in slope 

of the pairs of lines drawn through Cp* in figure 12. 


The experimental results also indicate that 


d 2 Mg 

dMj 


M^l.OO 


= 0 
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This then implies (see appendix B) that 

d2c P _ 127 + 4 | lOy + 6 c ^ 

dM °° 2 Mo<,=i.oo + x ) 2 (7 + 1) 2 D 

and an estimation can then he made of the range about = 1.00 where 
equation (l) is valid. 


Shock -Wave Detachment 

Conical -flow theory indicates that for a given cone angle of a 
semi -infinite cone a certain minimum Me is reached below which a coni- 
cal solution is no longer possible. This value of is defined to be 

the shock-wave detachment . However, a finite cone introduces a char- 
acteristic length into the problem so that curved attached shock waves, 
which would provide the necessary pressure gradient to turn the flow near 
the cone tip, may exist at values of less than the conical detach- 

ment Moo • 


Present experimental results indicate only that shock-wave detach- 
ment for a given cone angle does not occur at an greater than that 

predicted by conical theory. A collection of data from references 22, 

25, 26, and 15 is shown in figure 13. The ratio s/d where 6 is the 
center-line distance from the shock wave to the plane of the cone shoulder 
and d is the body diameter at the shoulder, that is, the sonic point, 
is seen to approach asymptotically to the value of 6/d at attachment. 

The asymptotic behavior of 6/ d complicates the fairing of the proper 
curve of 6/d versus M^ particularly in view of the paucity of experi- 
mental points in the immediate vicinity of shock-wave attachment. 

inference 27 data are presented for the shock-wave detachment 
distance of several cone angles at M^, = 2.45* A discrepancy was found 
between the experimental and theoretical values of the cone angle at 
which shock-wave detachment occurs for a fixed value of M^, detachment 
appearing to occur at a cone angle slightly greater than that predicted 
by the conical theory. This behavior would correspond to shock -wave 
detachment for a fixed cone angle occurring at a value of less than 

the theoretical conical value of M^. Again, however, the discrepancy 
may be caused by the manner in which the experimental curve was faired. 

Thus , the experimental results appear to indicate only that shock- 
wave detachment for a finite cone occurs at a value of either less 

than or equal to the value of M^ predicted by conical theory but not 
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at a larger Mco. The fact that shock -wave detachment does not appear 
to occur at a value of M^, greater than that predicted by the conical 
theory indicates that the presence of a boundary layer on the cone tip 
does not affect the conicity of the flow near the cone tip to the extent 
of precipitating detachment of the shock wave. 

Figure 13 also indicates that when the shock-wave detachment dis- 
tance is large, the position of the shock wave is dependent only on the 
diameter of the cone at the shoulder or sonic point and is independent 
of the cone angle. When the shock wave is quite near the cone tip, how- 
ever, the detachment distance is also dependent on the cone angle. This 
manner of behavior of the shock-wave separation distance was discussed 
by Busemann (ref. 28) and was shown experimentally for two-dimensional 
wedge sections by Griffith (ref. 29)- 


Transonic Similarity 

The transonic -similarity rules for the drag coefficient and pres- 
sure coefficient, as derived in reference 10, cannot be checked by the 
experimental results of this report. The derivation assumes that the 
cone -surface boundary condition is the approximate tangency condition 
which is valid for relatively small angles . A 20° semiangle cone is 
the minimum -angle cone for which detached shock-wave flow can be obtained 
in the transonic wind tunnel, and thus the experimental models were 20°, 
25°, and 30° semiangle cones. The experimental cone angles are much 
larger than the cone angles for which the approximate tangency condition 
is reasonable, and, therefore, the transonic -similarity rules of refer- 
ence 10 are not applicable . 


THEORETICAL CONSIDERATIONS 


At the present time, theoretical solutions have not been found to 
describe the flow about a finite cone for the complete Mach number range. 
Theoretical solutions are available for only two Mach number regimes, 
namely: 

(1) Exact conical theory may be applied if the surface Mach number 

is greater than sonic 

(2) At Mca = 1.00, Yoshihara (ref. 11) has calculated the flow about 

small -angle cones by a relaxation technique 

No solution has been determined if Moo is less than sonic. How- 
ever, Van Dyke's second-order supersonic theory and technique (ref. 30) 
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possibly can be applied to the subsonic case since, if the appropriate 
changes of sign are made in the particular solution found by Van Dyke 
for the supersonic case, the particular solution becomes valid for the 
subsonic case. 

A solution remains to be found for the regime between Moo = 1.00 
and the value of M^, where M g becomes equal to 1.00. The problem 

would be greatly simplified if the transonic equations could be employed. 
To test the feasibility of the approximations inherent in the transonic 
equations, an approximate solution has been found for conical flow about 
a semi-infinite cone using the transonic equations. The details of the 
solution are presented in appendix C . 

The solution is compared with the exact conical theory in figures 14 
and 15 . Figure 1^ shows the comparison between the shock-wave angles 
predicted by the exact theory and by the transonic approximation. The 
surface Mach number comparison is shown in figure 15 . From figure 15 it 
can be seen that the transonic approximation is quite satisfactory and 
is probably better than slender-body cone theory, since slender-body cone 
theory does not consider the presence of the conical shock wave. Also, 
figures Ik and 15 show the agreement of the cone single at shock-wave 
detachment as found from the transonic solution and from the exact theory. 

The above comparison of the exact conical solution and the approxi- 
mate solution indicates that the transonic equations contain all the 
terms of importance in the exact equations for conical flow about cones, 
so that the transonic equations may be employed with confidence in the 
range of M^ from M*, = 1.00 to an M^ for which M s = 1.00. 


SUMMARY OF RESULTS 


The following results were obtained from an investigation of tran- 
sonic flow past cone cylinders . 

1. The variation of drag coefficient Cp with free -stream Mach 
number M^ was determined experimentally. The slope of Cp versus Mo, 
at Moo = 1.00 agrees with the theoretical prediction. The deviation 

of Cp versus Moo from the conical flow value of Cp when M s < 1.00 
is demonstrated. 

2. The experimental results for the shock -wave angle, sonic-line 
location, and surface Mach number in the region near the cone tip indi- 
cate that the flow is conical near the tip of a finite cone even when 
the surface Mach number is less than sonic. The surface Mach numbers 
for the rest of the cone deviate from the exact conical values when 
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Ms < 1.00. Also, a case of shock-free supersonic -to-subsonic compression 
is demonstrated experimentally. 

3. An approximate solution for transonic conical, flow has been 
developed and the agreement with the exact conical theory indicates that 
the axially symmetric transonic equations retain the important features 
of the exact equations. 

A. Present experimental values of the detachment distance of a 
shock wave from a finite cone tip do not demonstrate agreement with the 
detachment Mach number predicted by conical theory for a semi -infinite 
cone and the question of shock-wave detachment from a finite cone remains 
undecided. 


California Institute of Technology, 

Pasadena, Calif., April 15, 1953- 
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APPENDIX A . 

REDUCTION OF INTERFEROMETER DATA 


An interferometer determines the advancement or retardation of a 
light wave in a medium with respect to a coherent light wave in a ref- 
erence medium. Since the wave velocities are a function of the indices 
of refraction of the respective mediums and consequently of the densi- 
ties of the mediums, it may be shown that 



(Al) 


where 

p (l) density of undetermined medium 

density of reference medium 

l path length 

k Glads tone -Dale constant 

wave length in vacuum of light employed 

S fringe shift; in case of finite -fringe interferograms this 

is ratio of displacement of a fringe to interval between 

undisturbed fringes 

In equation (Al) it has been assumed that the light beams traverse 
identical geometrical paths, so that refraction, if present, is neglected. 
Also, the relationship between the index of refraction n and density is 
assumed to be 


n = 1 + Kp 


(A2) 


If n = 1 + a where a « 1, equation (A2) is obtained by linear- 
izing the Lorentz law. 


For the axially symmetric case, the fringe shift for a light path 
perpendicular to the axis of symmetry becomes 




r dr 


(r*-y2) 1/2 
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where r is the radial distance from the axis of symmetry to a point on 
the light path and y is the perpendicular distance from the axis of 
symmetry to the light path. 

In the present investigation, the density field was bounded by a 
shock wave at a distance y g from the axis and the reference density 
was the free -stream density p ; thus, 

OO 7 


S(y) . &L Ifoj - pjr dx (AJ) 

( r 2 . y 2) l/2 

Weyl (ref. 31) introduced the assumption that S(y) is a linear func- 
tion of in a small interval of y. The validity of this assumption 

for the present investigation is indicated by the parabolic nature of 
the typical fringe-shift curves shown in figure 16. If the substitutions 


v = v'- 


u 




u s = y c 


are made, equation (A3) becomes 


k r u s [p(v) - pjav 


s (u) = f f 

A o J u 


l/v^“u 


(A4) 


This is analogous to the solution of Abel's problem. The solution of 
equation (A4) for p is 


P(w) " Poo = - 

00 . rt 


r Us (dS/du)du 
' w |/u - W 


(A5) 


where w = r . A complete proof of the solution may be found in refer- 
ence 19. 
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If the region of integration in equation (A 5 ) is divided into equal 
increments in r of width b. 


r 


i 


ib 


where i is an integer. A numerical evaluation of equation (A 5 ) is then 


rtKa Pco k=i 


(a6) 


where 


r s = Nb 


and 


i/(k + l) 2 - i 2 - l/k 2 - i 2 


The above solution of the problem is essentially that of reference 18 
and a table of A^q for 50 intervals will be found in reference 18. 

From the density ratio determined by equation (A 6) the local Mach 
number may be computed. An approximate correction to the local stagna- 
tion pressure p Q ' downstream of the shock wave was made by assuming 

that on the cone surface p Q ' was determined by the nose shock-wave 
angle and at a given chordwise station p D ' varied linearly with the 

value on the shock wave at the given station. If the approximate stream- 
lines are calculated, as in reference 20, a more refined correction is 
obtained. From the experimental values of the local Mach number a topo- 
graphic map was plotted, and from this map the desired Mach number con- 
tours are found. 


Conical Fields 

If it is desired to determine whether or not a given axially sym- 
metric field is conical, a simple test can be made. Returning to equa- 
tion (A3) it is assumed that 


P = p(n) 
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where q = r/x and x is the axial distance from the conical origin. 
Then equation (A3) becomes 


s( y ) = f 

x A / 



% [?(*)) - pJ 1 ! dT l 

,1 ^2 



where tj is the tangent of the shock -wave angle and tl = y/x. Thus, 
s 1 

if the field is truly conical. 


s(y) = f 



X 


and a plot of S(y)/x versus y/x for various values of x will yield 
a group of coincident curves. Examples of this technique are shown in 
references 19 and 20. It is interesting to note that in reference 20 
figure 6 (b) indicates conical flow near a cone tip for flow regime IV, 
that is, a 35 ° semiangle cone at M*, = I. 87 , using the above technique. 


Model size .- From equation (A3) it is evident that the fringe shift 
at a particular chordwise and radial station is a linear function of the 
model size for fixed values of density. Thus the model should be as 
large as is compatible with the test -section dimensions with regard to 
blocking and so forth. 

Finite fringe spacing .- The fringe spacing in the undisturbed field 
must be such that a sufficient number of data points may be determined 
between the shock wave and the cone surface at a particular chordwise 
station. However, for a given fringe shift S the displacement of the 
fringe is proportional to the undisturbed fringe spacing, and the accur- 
acy of the fringe data will be improved by increasing the undisturbed 
spacing. A compromise must be effected between the desire for many 
fringe shift points at a given chordwise station and the accuracy of the 
individual points. In the present investigation, this compromise pre- 
cluded a study by interferometry of the flow properties in the immediate 
vicinity of the cone tip when the shock wave was attached. 

Accuracy . - The accuracy of the interferometric method is affected 
by refraction, inhomogeneities in the reference flow, the numerical 
approximation, and so forth. An estimate of the accuracy can be obtained 
by noting the comparison of interferometer data and shock-reflection data 


Subsidiary Considerations 
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in figure 3 and the values of local Mach number behind the shock wave 
indicated in figures 6 and 7. 

Comprehensive discussions of the accuracy of the method may be 
found in references 18 and 19. 
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APPENDIX B 

VARIATION OF PRESSURE AND DRAG COEFFICIENTS NEAR = 1.00 


The pressure coefficient 


C 


P 


2 / Ps 

rHx > 2 \ p °° 


where p is the surface static pressure and. p is the free-stream 

S co 

static pressure may, in transonic flow, be written as 


C 


P 



- 1 


(Bl) 


where M s is the surface Mach number. In equation (Bl) it has been 
assumed that 




that is, the stagnation pressure loss across any shock waves may be 
neglected. 

Defining £ and | as 


i = M «. 2 


1 



Equation (Bl) becomes 


c P = f(e,s) 


2_ 



1 + 


7 - 1 


V 



(B2) 
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For a fixed body geometry, M s is a function of only, 

therefore 


I - !(?) 


Thus, the first and second total derivatives of C p with respect to £ 
are 



f f + 


1 dC 


(B 3 ) 


and 




+ 2f 


51 







(B4) 


In general, d|/d£ and d^|/d£ 2 are not known , however, an argument 
presented by Liepmann and Bryson (ref. 24) shows that, for {; = 1.00, 
d|/d^ = 0. The same argument, namely, that M s has a stationary value 

at Moo = 1.00, cannot be used to evaluate d^|y^d(^. However, an inspec- 
tion of the experimental data (see figs. 10 and 11) indicates that the 
curve of M s versus M^ has an inflection point at Moo = 1-00 and, 

thus, d 2 £/df; 2 = 0 at t; = 1.00. 


Using equations (B 3 ) and (B4) and the above argument, the deriva- 
tives of c p become 



where ( )* indicates evaluation at 5= 1.00, that is, at Moo = 1.00. 
The derivatives of C p with respect to Moo are then 



2f 



= 4 f££* + 2f 
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The partial derivatives f^ and f may be calculated from equa- 
tion (B2) and thus 


Fv) 

IdMooi 


7+1 7 + 1 


(B 5 ) 


\dMocw 


_ 127 + 4 IO 7 + 6 * 

(7 + l ) 2 (7 + l ) 2 P 


(b6) 


The fore drag coefficient of a finite cone is given by 


C D 



where a = ^r/r-^j , r is the local cone radius, and r^ is the cone 
base radius. Using equations (B 5 ) and (B6), the derivatives of Cp are 


and 



4 

7+1 


2 

7 + 1 


C D * 


d 2 C D \ _ _ 127 + 4 + 107 + 6 c # 

^dMoo 2 / (7 + i ) 2 (r + i ) 2 


(B 7 ) 


(B8) 


A comparison of equation (B7) and experimental data is shown in 
figure 12. Equations (B5 ) and (B 7 ) were given previously by Bryson 
(ref. 4). 

It should be noted that the first term 4/(7 + l) of equation (B5) 
may be derived from the linearized transonic Cp (ref. 1) which is 


C 


P 


= 2 


K> c 



(7 + 1 ) 


(B 9 ) 
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The second term -[ 2/(7 + l)jCp* of equation (B5) is then of the nature 
of a second-order correction term. However, if (d 2 C p /dMoo 2 )* is com- 
puted from equation (B9), the result does not agree with the first term 
of equation (b 6 ) and is, in fact, of opposite sign. 
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APPENDIX C 

TRANSONIC APPROXIMATION FOR CONICAL FLOW 


In discussing axially symmetric transonic flow, the following 
approximations to the exact equations and boundary conditions are employed. 

If q x is the velocity in the axial or X-direction and q r is the 
velocity in the radial or r-direction it is assumed that 


a* + 0 x 



where a* is the velocity of sound at M^ = 1.00, 0 is the perturba- 

tion potential, and $ x , 0 r « a*. Then, defining 


u 


7 + 1 
a* 


K 


v 


7+1 

a* 


0r 


the continuity equation is approximated by 


uu x - v r " p v = 0 ( C1 ) 

and the condition of irrotationality becomes 

u r - v x = 0 (C2) 

The exact tangency condition on the body surface is replaced by 

v = (7 + 1) tan 0 
o ' ' 


(C3) 
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where v Q is v on the surface of the body and 0 is the inclination 
angle of the body surface . The above relations are derived in greater 
detail in reference 9 . 

To test the feasibility of the approximations inherent in equa- 
tions (Cl), (C2 ) , and (C3) an approximate solution for the flow about 
a semi -infinite cone will be developed. This approximation will then 
be compared with the Busemann solution of the exact equations . 


Conical Solution 


Assuming that u and v are functions of a where 


r 


equations (Cl) and (C2) become 


uu a + av a - v = 0 


(Cla) 


au 0 + v a = 0 


(C2a) 


A solution of the form 


u = f(v) 


(C4) 


will be sought. From equation (C4) 


du _ du dv 
da dv da 


dv _ u du 
da da 


but from equation (C2a) 
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Therefore 


du _ _ 1 
dv a 

The relationship in equation (C 5 ) is exact (see ref. 32). 

Differentiating equation (C5) with respect to a yields 


dv = 1 

^- a p ^ u 

dv 2 

and therefore 

du _ _ 1 

^ -Z) d 2 u 

dv 2 

Substituting equations (C5), (C6), and (C7) into equation (Cla), 
result is 


d 2 u du _ 
dv 2 dv 



3 


The following conditions are imposed on equation (C8). 

(a) At the shock wave: If u g and v g are the values of 

at the shock wave, then, from the transonic -shock polar relation 


= (> - u s) \[ ] 


U«, + Uc 


where 


Uoo = Mco 2 - 1 


The shock -wave angle p is given by 


v 


s 


Uoo " u s 


(C5) 


(C 6) 


(C7) 

the 

(C8) 


u and v 

(C9) 


cot p = 
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hence, from equation (C 5 ), 




(CIO) 


(b) On the Cone Surface: The tangency condition must be satisfied, 

hence 


v Q = (7 + 1 ) tan 0 


(Cll) 


where v Q is v on the cone surface and 0 is the cone semiangle . 
From equation (C 5 ) 


= -tan 0 (C12) 

o 

An exact solution of equation (C8) has not been found. However, if the 
right-hand side of equation (C8) is assumed to be small, an iteration 
solution may be found. 


du 

dv 


First Approximation 

As a first approximation to the solution, set the right-hand side 
of equation (C8) equal to zero; that is. 



The solution of equation (C13) is 


(C13) 


u = A log e Bv 


(C14) 


Applying conditions (C9) and (CIO) to this solution, equation (Cl4) 
becomes 


U = U 3 - K - U s) l0 S e ^ 


(C15) 
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Second Approximation 

As a second approximation to the solution of equation (C 8 ), the 
right-hand side of equation (C 8 ) is assumed to be given with sufficient 
accuracy by equation (C15)- Hence 



(C16) 


The solution of equation (Cl 6 ) is 


u = 


_ (uco - u s ) 3 

^V 2 


" u s + - u s) + K - u s) lo S e 7 - 

s 


+ C log e Dv 

(C17) 


Applying conditions (C9) and (CIO) to this solution, equation (C17) 
becomes 


u = K - u s) 3 
Av 2 


-U S + (Uoo - U S ) + (Ut» - Ug) !Og e JE- 


(Upo + 3 u s) - Ug) log _v_ _ Upo(“po - 5u s ) 


2 (tloo + Ug) 


V s 2 (uoo + u s ) 


(C18) 


The values of u s and v s appearing in equation (Cl 8 ) are not 
independent of the cone semiangle 0. Applying condition (C12) to 
equation (Cl 8 ) and solving for v s , the result is 


l°g e v s 


= l°g e (7 + !) tan 9 + 



+ 1) tan 9 ] 2 (uqq + 9u s ) 
(Uco + u s) K. - %) 5 



2(7 + 1)^ tan^9 


K - u s ) 4 


(C19) 
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The transonic -shock polar relation also must be satisfied, hence 



+ u s 


(C20) 


Because of the nature of equation (C19), explicit solutions of equa- 
tions (C19) and (C 20 ) for u s and v s in terms of u^, and 0 have 
not been found. However, a solution may be found graphically. 

The values of u s and v s thus determined for a given value of 0 
and Uo, may be introduced into equation (Cl 8 ). If the value of u on 
the surface u Q is desired, then substitution of 


in equation (Cl 8 ) yields u Q . 

It should be noted that, for a given value of 0 and u^,, two sets 
of values for u s and v s are found. These correspond to the "strong" 
and "weak" shock waves . 

Comparison of Second Approximation and Busemann Conical Solution 

Wave angle . - The wave angle p is determined by the values of u„ 

s 

and v s since 


The degree of agreement is apparent in figure 14. 

Surface Mach number .- In terms of u Q and v 0 , the surface Mach 


V 


v 0 = (7 + 1 ) tan 0 


number is 



(C21) 


where 
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The usual transonic approximation reduces equation (C21) to 


M s 2 = 1 + uq 


(C22) 


equation (C21) is shown in figure 15 • The agreement with the exact 
theory is quite satisfactory. 

Apple curve .- Two sets of values of u g and v s will satisfy equa- 
tions (C19) and (C20). The two solutions correspond to the "strong" and 
"weak" shock-wave solutions predicted by Busemann. If the values of u Q 
and v Q for all possible cone angles and for both types of solution are 
plotted in the hodograph, the resultant curve is Busemann' s "apple" curve. 
The apple curve found in the second approximation is shown in figure 17. 
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Figure 1.- Finite-fringe interferogram for 30° 

^ = 1.280. 


L-8I1.893 

semiangle cone cylinder. 
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Figure 2.- Wave reflection from sonic line for 20° semiangle cone cylinder. 

^ = 1.297. 
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Figure 3-- Sonic-line location for 25° semiangle cone cylinder. 


NACA TN 3213 


39 



Figure 4. 


- 20° semiangle cone cylinder , 


L-8I4.895 
Me = 0.942. 
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(a) Moo = 1.328. 

Figure 5*- 25° semiangle cone cylinder. 
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(b) Moo = 1.415. 
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Figure 6.- Mach number field for 25° semiangle cone cylinder. 
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(b) Mco = 1.253. 
Figure 6.- Continued. 
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M m = 1.273 



r / r b 


0 .5 ~ L0 


(c) Moo = 1.273. 
Figure 6.- Continued. 
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Figure 6.- Concluded. 
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Figure 8.- Sonic-line location by shock reflection for 20° semiangle 

cone cylinder. 
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Figure 10.- Surface Mach number on 25° semiangle cone cylinder. 
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Figure 12.- Drag coefficients for cone cylinders. 
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Figure 13 .- Shock- wave detachment distance for finite cones. Data are 

from references 15, 22, 25, and 26. 
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Figure l4.~ Shock-wave angle for semi- infinite 


cones, = 0.25; M m = 1.12. 



Figure 15 •- Transonic conical solution; surface Mach number versus cone 
semiangle. Uoo = 0 . 25 ; Moo =■ 1.12. (Moo is defined as Moo = s/T~T~uZ.) 
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Figure 16.- Interferogram fringe shifts for 25° semiangle cone cylinder. 
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Figure 17.- Transonic apple curve at u^ = 0.25 from second approximation. 
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